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Recent experimental results suggest that the neutrinos of the Standard Model are massive, though light. 
Therefore they may mix with each other giving rise to lepton flavour or even lepton number violating processes, 
depending on whether they are Dirac or Majorana particles. Furthermore, the lightness of the observed neutrinos 
may be explained by the existence of heavy ones, whose effects on LFV would be very sizeable. We present an 
analysis of the effect of massive neutrinos on flavour-changing decays of the Z boson into leptons, at the one-loop 
level, independent of neutrino mass models. Constraints from present experiments are taken into account. 



1. Introduction and motivation 

With the Giga-Z option of the Tesla linear 
collider project one may expect the production 
of about 10 9 Z bosons at resonance [ EJ. This 
huge rate, about a factor 100 higher than at 
LEP 1, allows one to study a number of prob- 
lems with unprecedented precision. Among them 
is the search for lepton-flavour changes in Z de- 
cays: Z — > e/i, /it, er. The best direct limits are 
obtained by searches at LEP 1 (95% c.l.) [f|: 
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A careful analysis shows that the sensitivities 
could be largely improved at the Giga-Z [ ||] , 
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Non-zero rates are expected if neutrinos are 
massive and mix [ ^| . From experiments we have 
evidence of tiny neutrino masses and substantial 
mixings. Unnaturally small mass scales may be 
indicative for a mechanism which produces at the 
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same time very large masses. Heavy neutrinos 
are introduced in most GUTs [ || and string- 
inspired models [ |j, and are suggested by the 
seesaw mechanism [ ^| . 

The above observations motivate us to have a 
closer look at the prospects of observing LFV. We 
will explore here the following scenarios: 

(i) The z^SM. We treat the known n G = 3 
generations of light neutrinos (v e , v^, v T ) as mas- 
sive Dirac particles. Individual lepton numbers 
L e ,L^,L T are not conserved (in analogy to the 
quark sector). As a by-product, the Z decay am- 
plitude into two quarks of different flavours can 
be read off from our general expressions. 

(ii) The z^SM extended with one heavy ordinary 
Dirac neutrino (usual SU(2)x,®U(l)y quantum 
numbers). This case implies the existence of a 
heavy charged lepton as well. It is a simple ap- 
plication of case (i) for heavier neutrinos. Again, 
total lepton number L is conserved. 

(iii) The i/SM. extended with tir = 2 heavy 
right-handed singlet Majorana neutrinos. Not 
only individual, but also total L is, in general, not 
conserved since the presence of Majorana mass 
terms involves mixing of neutrinos and antineutri- 
nos, with opposite lepton number. For two equal 
and heavy masses this case reduces to the addi- 
tion of one heavy singlet Dirac neutrino. In this 
latter case L is recovered [ |§| . 

Our results are independent of neutrino mass 
models. We take into account constraints on neu- 
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trino masses and mixings given by oscillation ex- 
periments (light sector) or imposed by unitarity 
and precision tests of flavour diagonal and nondi- 
agonal processes (heavy sector). 

2. Z — > with massive neutrinos 

The amplitude for the decay of a Z boson into 
two charged leptons with different flavour (l\ and 
£2) vanishes in Born approximation but receives 
quantum corrections. It can be written for mass- 
less external fermions as: 

M = -^^V £ ^ 2 (p 2 )7^1-75)^(-Pi) (7) 

with aw = a/s w . The dimensionless form factor 
V is a function of Aq = Q 2 /M^, with Q 2 = (p 2 — 
Pi) 2 (to be fixed at the Z peak), and the masses 
Xi = m^/My^ and mixings of the massive, virtual 
neutrinos i [ [| [l0|. V gets contributions from 
one-loop vertex- and self-energy-graphs, that are 
calculated here in the 't Hooft-Feynman gauge. 
The former consist of triangle graphs: two virtual 
neutrinos coupled to the external Z boson with a 
W or a Goldstone boson <f> being exchanged (W, 
</>); and one virtual neutrino exchanged with W 
or cj> coupled to the Z (WW, <j><j> and Wcf>). The 
self-energy graphs are corrections to the external 
fermion legs (£). 

One must discuss separately the cases of ordi- 
nary Dirac and general Majorana neutrinos [ |TT): 

Dirac: 



(8) 



V(i) = v w {i) + v^i) + v ww (i) + v^i) 

+vw4>( i ) + v-z{i), (9) 



Majorana: 



n G +n R 



V M = E B ^ B tjVM(i,j), (10) 

VAiihj) = v w (i,j) + v<j,(i,j) + vww{i) 

+v < j >< j ) {i) + v w <f>(i) + v-s(i), (11) 

where only global mixing factors are extracted: 
the leptonic CKM mixing matrix V for ordinary 
Dirac neutrinos, and its generalized version, B, 



for Majorana neutrinos. The latter mixings ap- 
pear in charged-current interactions of the left- 
handed (LH) components of uq generations of or- 
dinary charged leptons (f L = cl, ■ ■ •) with 
n<3 LH isodoublet neutrinos (v° L = v e , v^, v T , . . .). 
The interaction eigenstates are in general not the 
same as the physical mass eigenstates but a mix- 
ture [HI HII, 

no riG+n R 

^ = E U &^> v li= E v ijVLj , (12) 

j'=i i=i 
where v—n v c are uq + Ur Majorana fields (self- 
conjugate up to a phase rf). In the physical basis, 
9_ 
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C cc = -^W^ Li ^P L u Li + h.c. 
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where P R . L = \ (l ± 75) and 
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is a rectangular uq x (uq x ur) matrix. 

The main feature distinguishing Dirac and Ma- 
jorana cases is the existence of nondiagonal ZviVj 
vertices (flavour-changing neutral currents), cou- 
pling both left- and right-handed components of 
the Majorana neutrinos to the Z boson, 
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a quadratic (no + «_r) 2 matrix. Such vertices 
appear in graphs where a W or a Goldstone boson 
<f> is exchanged: 

vw(i,j) = -C y - [Aq(C + Cn 1 
-2C 24 + 1] 

+ C* - y/XiXj 
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The (diagonal) contributions of these graphs for 
Dirac virtual neutrinos [vw(i) and v^,{i,j)\ are 
obtained by the replacements: 



c 



► C..(Xi,Xi), (19) 

ij — * {vi + Oi) Sij = Sij , (20) 

C£ -> -(«,•- o<) =0. (21) 
The rest of the contributions are: 

v ww (i) = 2c^ (27^ ) [A Q (Cn + C« + C 23 ) 

-6C 24 + 1], (22) 

v^(i) = -(l-2s 2 w ) (27^) Xi C 24 , (23) 

«wr^(i) = -24, (27^) Aj Co, (24) 

Vs{i) = -(«i + fli-44rOj) 

x [(2 + Ai)Bi + l]. (25) 

We have introduced above dimensionless two- 
and three-point one-loop functions: 

Bj = M 1 (X i ) = B 1 (0;m 2 i ,M^), (26) 
C. = C..(Ai) 

= M^C.XQ,Q 3 ,0',mlM^,M^), (27) 
C. = C.XXuXj) 

= M 2 W C. (0,Q 2 ,0;M^, m 2 ,™ 2 ), (28) 
from the usual tensor integrals [ |l5) 

B"(p 2 ;mg,m?) = p"Bi, (29) 
C»(p 2 ,Q 2 ,p 2 ;m 2 , m 2 ,m 2 ) = P %C n +p%C 12 , (30) 
C» v (pi,Q 2 , P l-,mlm 2 ,m 2 ) 

= P >tplC 21 + p^ 2 C 22 + +^PDC23 

+ 5 ^M^C 24 . (31) 

The tensor integrals are numerically evaluated 
with the computer program LoopTools [ 

The expressions above, in the Dirac case, are 
also valid for quark flavour-changing Z decays in- 
volving virtual quarks i, just by substituting their 
corresponding weak isospin , electric charge , 
and vector u$ and axial- vector a, couplings, 



v>i — 1^ 2eiSyy, (32) 
a, = 4 L - (33) 

It turns out convenient to cast (|TT| ) as 
V M (i,j) = SijFiXJ + CijGiXuXj) 

+C* ij ^X~X j H(X i ,X j ). (34) 



The Dirac form factor (0) is then 

^(i)= J F(Ai) + G'(A i ,A < ). (35) 

The amplitude .M is finite without renormaliza- 
tion, but the form factors V and Vm are not. The 
divergences are such that they exactly cancel due 
to unitarity relations among the mixing matrix 
elements V, B and C[||,0. 
The branching ratio Z — > ifl^ reads 



a 3 w M z 



(36) 



The constant in front of pj) is of O(10" 6 ) 



3. Contribution from light neutrinos: i/SM 

Ignoring the not yet confirmed results of the 
LSND accelerator experiment, all neutrino ex- 
periments are compatible with the oscillation be- 
tween two of a total of three neutrino species. The 
signals in atmospheric experiments [ |Ti| are at the 
90% c.l. compatible with — v T oscillations, 
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Solar experiments [ Il8| indicate v e — mixing, 
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10~ iu -h 10 _o eV z , 
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(39) 
(40) 



(there are solutions for vacuum and matter os- 
cillations compatible with a wide range of mixing 
angles). From reactor searches, there are no hints 
of v e — v T oscillations [ |l9) , which implies 



sin 2 26»i3 = 0. 



(41) 



Taking this information into the standard param- 
eterization for the mixing matrix [ |^| (oscilla- 
tion experiments are insensitive to Majorana CP- 
phases) and putting the Dirac CP-phase 5 = 0, 
since no information on it is yet available, one has 
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Using the unitarity of V, with l\ ^= £ 2 , 

v3 



BR(Z -> efif) = 
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192^c 2 w T z 

3 

E v « v u*w-^(°)] 



(43) 



Performing a low neutrino mass expansion of the 
tensor integrals (A, -C 1) one finds [|2(], |2l) 

V(Xi)-V(0) = ai X t + 0(X 2 ), (44) 
ai = 2.5623- 2.950 i. (45) 

Therefore BR(Z -> ifef) goes as mf for low neu- 
trino masses. This behaviour is still a good ap- 
proximation not far below the Z mass (see Fig. 1). 

Substituting the phcnomcnological squared 
mass differences A« = Am|/M^, and the mixing 
angles, one has [ O 
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4 x 10~ 60 , (46) 
3x 10" 6 x | S ? 2 A 12 -A 23 | 2 
(3 -=-30) x 10" 55 . (47) 



BR(Z — > /i T r ± ) 



These rates are extremely small. In fact, the con- 
tributions from the observed light neutrinos will 
be neglected in the next section, where we extend 
the vSM to accommodate heavy neutrinos, taking 
massless the light ones. 

4. Contribution from heavy neutrinos 

Assume that there is a sector of heavy neutri- 
nos Ni mixing with the light ones. Using a general 
formalism developed in [ ^2] one can exploit mea- 
surements of flavour diagonal processes (checks of 
lepton universality and CKM unitarity, Z boson 
invisible width, ...) [ |23| to obtain indirect exper- 
imental bounds on light-heavy mixings [ |l4| : 



(48) 



(replacing B by V for Dirac neutrinos) . The most 
recent indirect bounds [ |24j] : 

s 2 Ve < 0.012, si < 0.0096, s 2 Vt < 0.016, (49) 
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Figure 1. Contribution of one massive Dirac neu- 
trino (solid) or up-quark (dashed) to BR normal- 
ized to unit mixing. Dotted lines using Q 2 = 0. 



are only improved by direct searches for flavour 
nondiagonal processes involving the first two lep- 
ton generations. In fact, from BR(/i — ► cy) < 
1.2 x 10" 11 [H| one may infer [|§]: 

< 1.4 x 10~ 8 . (50) 
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4.1. i/SM + one heavy ordinary Dirac 

From unitarity of V for uq + 1 generations, 



BR{Z 



l92ir 2 c 2 T z 



x \V(X N )-V(Q)\ 2 . (51) 



The form factor V, subtracted, squared and nor- 
malized is depicted in Fig. 1. The results agree 
with earlier calculations [ also for the quark 
flavour-changing Z decays [|10|. As expected, the 
approximation Q 2 = is very bad for tojv < Mz 
but it makes sense for mjv > Mz (Fig. 1). 
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/i T r ± in the i/SM extended with: (i) one 



Figure 2 
of Z 

heavy ordinary (thick solid) or singlet (thin solid) 
Dirac neutrino of mass m n 1 ; (ii) two heavy right- 
handed singlet Majorana neutrinos (dashed lines) 
with masses mjvi and m^ 2 . The upper limits of 
19|) are taken as light-heavy mixings. 



4.2. i/SM + (nji — 2) Majorana neutrinos 

Unitarity constraints on B and C [ [l4| allow to 
write Vm m terms of the heavy sector only: 



V M = 



F(A JVi )- J F(0) + G(A^,0) 

+<7(0,Ajv 4 ) - 2G(0,0)] 
+0^^ [ G(\Ni , A at,. ) — G(\Ni , 0) 
-G(0,A^) + G(0,0)] 

+Cjv' i jv i y^Ni^Nj H(X Ni , X Nj )^. (53) 

For n_R = 2 the mixing matrices are exactly 
calculable in terms of and r = m 2 N ^/m 2 Ni [ 
fL4| . The upper values for the branching ratios 
can be then straightforwardly obtained from the 
bounds (p9|), given the heavy masses m^, 
(Fig. 2). The case itin 1 — mjv 2 is equivalent 
to one heavy singlet Dirac neutrino (in fact, two 
equal mass Majorana neutrinos with opposite CP 
parities form a Dirac neutrino). 

In the large neutrino mass limit (Ajv x 3> 1) one 
obtains [ Ell 



Vm(Ajvi)?"; s vt ) 
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A Q InAjv, +0(1). (54) 



£9|), Fig. 2 shows the 



Taking the present upper bounds of the mixing 
matrix elements from ([T~ 
maximal BR(Z — ► /i^r^), for illustration. 

It is worth noticing that the expansion of tensor 
integrals in the large mass limit [ ^ |ll| at the Z 
peak yields 



1, 



^(Aat)-^(O) = -[A w + 2 



In X N 

(6.99 + 2.11 i)] + 0(lnX N /X N ), (52) 



leading again to an m 4 ^ growth of the branching 
ratios for large neutrino masses. 



The constant in front of the In A term is identical 
to the one in the Dirac case (|52]). In the approx- 
imation Xq = Q 2 /M^ — 0, the expression Q) 
is in agreement with [ M , but we take the ac- 
tual value Aq = M§/M^T= \jc 2 w « 1.286 in the 
whole calculation. 

Finally, notice that large neutrino masses are 
restricted by the perturbative unitarity condition 
on the decay width of heavy neutrinos [ |27| , 

r N ^(2x)^-m 3 m J2\BiNf<lm Nt (55) 

for Dirac (Majorana) neutrinos. In this way, the 
unacceptable large-mass behaviour of the ampli- 
tudes (cx m 2 N ) is actually cured when a sensible 
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light-heavy mixing (at most oc m^f) is taken into 
account [ |2^]. The restrictions (gq) lead to the 
end-points in the curves of Fig. 2. 

5. Summary and conclusions 

The perturbative unitarity limit on the decay 
width of heavy neutrinos effectively prevents their 
nondecoupling, thus ensuring smaller light-heavy 
mixings for increasing heavy masses. Given the 
present indirect upper bounds on light- heavy mix- 
ings of 0(1O~ 2 ), this already sets indirect upper 
limits on LFV Z decays of about BR < 10~ 6 with 
Majorana neutrinos and above LEP 1 reach for 
Dirac neutrinos. 

We have presented the full one-loop expecta- 
tions for the direct lepton flavour changing pro- 
cess Z — > ^1^2 with virtual Dirac or Majorana 
neutrinos. We conclude that: (i) the contribu- 
tions from the observed light neutrino sector are 
far from experimental verification (BR < 1CP 54 ); 
(ii) the Giga— Z mode of the future Tesla lin- 
ear collider, sensitive down to about BR~ 10~ 8 , 
might have a chance to produce such processes. 
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